Interaction of W(CO)6 with Si02 Surfaces - A Density Functional Study 



O 
<N 

■^-> 
o 

o 

<N 



Kaliappan Muthukumar, Ingo Opahle, Juan Shen, Harald O. Jeschke, and Roser Valenti 

Institut filr Theoretische Physik, Goethe- Universitdt Frankfurt, 

Max-von-Laue-Strafie 1, 60438 Frankfurt am Main, Germany 

(Dated: October 25, 2011) 

The interaction of tungsten hexacarbonyl W(C0)6 precursor molecules with Si02 substrates is 
investigated by means of density functional theory calculations with and without inclusion of long 
range van der Waals interactions. We consider two different surface models, a fully hydroxylated 
and a partially hydroxylated Si02 surface, corresponding to substrates under different experimental 
conditions. For the fully hydroxylated surface we observe only a weak interaction between the 
precursor molecule and the substrate with physisorption of W(C0)6. Inclusion of van der Waals 
corrections results in a stabilization of the molecules on this surface, but does not lead to significant 
changes in the chemical bonding. In contrast, we find a spontaneous dissociation of the precursor 
molecule on the partially hydroxylated Si02 surface where chemisorption of a W(C0)5 fragment is 
observed upon removal of one of the CO ligands from the precursor molecule. Irrespective of the 
hydroxylation, the precursor molecule prefers binding of more than one of its CO ligands. In the 
light of these results, implications for the initial growth stage of tungsten nano-deposits on Si02 in 
an electron beam induced deposition process are discussed. 

PACS numbers: 68.43.-h,68.43.Fg,71.15.Mb,71.15.Nc 
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INTRODUCTION 



The growth of patterned transition metal nano- 
deposits with tunable electronic properties has been in 
the focus of recent research towards nanoscale device ap- 
plications^ ^. One of the promising routes to achieve 
this goal is the preparation of nano-deposits from gas 
phase precursor molecules like tungsten hexacarbonyl 
(W(C0)6) or cyclopentadienyltrimethylplatinum(IV) by 
means of focused ion or electron beam induced depo- 
sition (EBID) or laser induced deposition techniques^lSl. 
While experimentally substantial progress has been made 
in the fabrication of transition metal nano-deposits with 
properties reaching from insulating to metallic or even 
superconducting behavioi^, a detailed understanding of 
the growth processes on a microscopic level is still miss- 
ing. This is due to a multitude of complicated processes 
involved in the growth of nano-deposits from the gas 
phase - like the adsorption or chemisorption of precursor 
molecules, heat transfer and the interaction with primary 
and secondary electrons in the EBID process as well as 
fragmentation and recombination of fragments - which 
do not allow a one step simulation of the growth process 
within reasonable time. 

W(C0)6 is one of the most prominent precursor ma- 
terials used in the EBID process and has been deposited 
on various substrates like Si02^ ^, W(llO)^, and Ni(lOO)^ 
by using a variety of techniques. The chemical and phys- 
ical properties of W(C0)6 as a molecule are well studied 
both experimentally and theoreticall}'^ ^^. However, to 
the best of our knowledge, there haven't been any reports 
of electronic structure calculations that treat the inter- 
action of W(C0)6 molecules with the Si02 substrate. 
Therefore, in this work, we focus on constructing a rea- 
sonable model for the Si02 substrate and study the inter- 
action of a single W(C0)6 molecule with the Si02 sub- 



strate as one important step towards the simulation of a 
realistic growth process in EBID. 

For the simulation of the W(C0)6 interaction with the 
Si02 substrate we employ two different surface models, 
a fully hydroxylated (FOH) and a partially (POH) hy- 
droxylated Si02 surface. The FOH - Si02 represents 
a realistic model for a substrate prepared under non- 
vacuum conditions in the absence of irradiation, while the 
POH - Si02 corresponds to a substrate under the influ- 
ence of an electron beam, where a partial removal of OH 
groups from the surface is expected due to the interac- 
tion with th e beam , similar to other surfaces such as Ti02 
and 7-Al20^^^^. The interaction between the precursor 
molecule and the substrate is then studied within density 
functional theory (DFT) using the generalized gradient 
approximation (GGA) with and without van der Waals 
(vdW) interactions. The impact of vdW interactions, 
which were found to be important for the interact ion of 
organic molecules on metal surfaces in recent studies^^EHl, 
is discussed for the present system. 



II. COMPUTATIONAL DETAILS 

DFT calculations for bulk Si02, isolated Si02 sub- 
strates and W(C0)6 precursor molecules, as well as 
the combined substrate/precursor molecule system were 
perfor med u sing the projector augmented wave (PAW) 
met ho (P^^^ as implemented in the Vienna Ab-initio 
Simulation Package (VASP-5.2.11)'^^H24|^ GG A in the 
parametrization of Perdew, Burke and ErnzerhoP^^was 
used as approximation for the exchange and correla- 
tion functional. In addition, dispersion corrections ^^ * ^^ * 
were used that simulate the long range vdW interactions, 
which are expected to be important for the adsorption of 
the W(C0)6 molecules on the FOH - Si02 surface. 
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FIG. 1: (Color online) (a) Structure of bulk /3-cristobalite 
Si02 and (b) side view of the slab geometry with a (111) 
surface used for this study. The corresponding unit cell is 
indicated. 



For bulk Si02, the ideal /3-cristobalite structure^ with 
cubic Fd3m symmetry (Fig. IT] (a)) and the experimental 
lattice constant a = 7.16 A was used. To simulate the 
substrate-precursor interaction, slabs consisting of 2 to 5 
layers of Si02 with a hexagonal Si02 (111) surface were 
generated from the experimental bulk structure. The 
minimal slab geometry to simulate the hexagonal Si02 
(111) surface has in plane lattice parameters correspond- 
ing to the primitive lattice constants a^ = a/^/2 = 5.06 
A of the face centered cubic cell. This minimal geom- 
etry does however not allow for a sufficient separation 
of adjacent images when the molecule is placed on the 
surface. Thus we used a 3 x 3 super cell geometry in the 
plane (in terms of primitive lattice constants) for all slab 
calculations. 

For the FOH - Si02 substrate, all surface silicon 
atoms (top and bottom) were terminated with OH groups 
(Fig.fl](b)). In the case of the POH - Si02 surface some 
of the OH groups were removed from this slab. We con- 
sidered three cases depending on the precursor molecule 
orientation with OH- vacancy concentrations of 11%, 22% 
and 33%. Both the FOH - Si02 and POH - Si02 sub- 
strates were then optimized. In a next step, we con- 
sidered the W(C0)6 molecule in similar orientations for 
both classes of substrates and the corresponding geome- 
tries were then reoptimized with fixed in plane lattice 
parameters but otherwise no symmetry constraints and 



the energies have been compared after relaxation. All 
calculations which involve the bonding of the precursor 
molecules to the substrate were carried out in a 4 layered 
3x3 supercell geometry in the plane perpendicular to 
the surface, corresponding to a low coverage of precur- 
sor molecules in the initial growth stage. The choice of 
the 3x3 supercell geometry (with up to 256 atoms for 
the substrate/precursor system) allows for a distance of 
more than 8 A between two precursor molecules (when 
measured between the two CO ligands of the adjacent 
molecules), thus providing an optimal compromise be- 
tween computational effort and minimization of the in- 
teraction between the precursor molecules. In all cases, 
the distance between slabs (on the c axis) was kept large 
enough (30 A that provides a distance of 24 A on the 
substrate/precursor system) to prevent significant inter- 
actions between the adjacent images. 

All calculations were performed in the scalar relativis- 
tic approximation. A kinetic energy cut-off of 400 eV was 
used and all ions were fully relaxed using a conjugate gra- 
dient scheme until the forces were less than 0.01 eV/A. In 
the geometry optimizations for the molecule and the sur- 
face models the Brillouin zone was sampled at the F point 
only. For the final comparison of energies and the den- 
sity of states (DOS), a k-mesh of 4 x 4 x 1 has been used 
for the slabs showing a convergence of the total energy 
to about 10~^ eV/f.u. For bulk Si02 the k-point sam- 
pling was performed with a 16x 16x 16 Monkhorst-Pack 
grid. For all DOS calculations a Gaussian broadening of 
20 meV was applied. Spin polarization was considered 
for all calculations. Different spin states for each clus- 
ter/adsorbate system were checked (i.e. in each case the 
two lowest possible spin states) and only the results of 
the ground state are reported below. 

In addition, the W(C0)6 molecule was optimized in 
the Oh symmetry in a cubic box of 30 A and the struc- 
tural parameters were used to compare with the surface 
models. We also a nalyzed the electronic orbitals with 
Turbomole 6.(PEIJ using the def-TZVP basis sets, the 
same GGA functional as for the bulk calculations as well 
as the RI (Resolution of Identity) approximation. In this 
case, an effective core potential has been used for W, 
which provides 12 active valence electrons. 

The adsorption energy was defined as AE = £^totai — 

^substrate ~ ^adsorbate^ whcrC £^totab ^substrate 7 ^^^ 

^adsorbate circ the energies of the combined system (ad- 
sorbate and cluster), of the cluster, and of the adsorbate 
molecule in the gas phase in a neutral state, respectively. 



III. RESULTS 

A. Electronic structure of the Si02 substrate 

Si02 exists in a variety of crystalline and amorphous 
modifications. In most experimental studies amorphous 
Si02 substrates are used for the growth process, which 
are however very demanding for the computational sim- 
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FIG. 2: (Color online) DOS of bulk Si02 (top) in comparison 
to the OH-terminated 4 layer slab (bottom). Shown are the 
respective total DOS (solid line) and the partial O 2p contri- 
bution for the bulk (dashed line). Results are shown for GGA 
calculations without vdW corrections. 




FIG. 3: (Color online) (a) optimized structure of W(C0)6 
and one of its triply degenerate (b) HOMO and (c) LUMO 
shown as a representative 



energy are dominated by O 2p states with some admix- 
ture of Si states (not shown). The electronic structure 
of the 4 layer slab shows a reasonable convergence with 
respect to the bulk (note that there are some smaller dif- 
ferences due to finite size and relaxation effects as well as 
due to different k-point sampling, which should however 
not affect the discussions in the subsequent sections). 



ulations. Therefore we use in our study /3-cristobalite, 
which resembles most closely amorphous Si02, as a rep- 
resentative structure for amorphous Si02. Experimen- 
tal studies have shown that /3-cristobalite and amor- 
phous Si02 have a very similar local structure and ex- 
hibit similar physical properties such as density and re- 
fractive index^^. Following Refs. [33H36l we use the ideal 
/3-cristobalite structure with Fd3m symmetry. 

It has been reported that the silica surface may contain 
segments of surfaces resembling both the (111) and (100) 
faces of cristobalite^^'^. The Bravais-Donnay-Friedel- 
Harker (BDFH) method provides an estimate of the rel- 
ative growth rates of the possible faces of each crys- 
tal structure and the resultant morphology, which shows 
that there exists only one dominant plane for Si02 which 
is (111)^^. Therefore, the surface of the substrate is sim- 
ulated by the (111) plane of bulk /3-cristobalite. 

For a slab size of 4 layers we have observed a con- 
vergence in structural and electronic properties. In the 
optimized structure of 4 layer slabs with no symmetry 
constraints the bond lengths of Si and O in the inner two 
layers are found to be in the range of 1.62 - 1.63 A match- 
ing well with the experimental values of 1.61 A^SI and 
previously reported theoretical values^^. Fig. |2] shows 
the DOS of bulk Si02 (upper panel) and the DOS of the 
FOH - Si02 4 layer slab (lower panel). The upper valence 
bands in the region between -10 and eV below the Fermi 



B. Electronic structure of W(CO)6 in the gas phase 

The structural and electronic properties of W(C0)6 
have been widely investigated using local basis sets over 
the last two decadeg^^ESl. To compare the electronic and 
energetic properties of W(C0)6 both in the gas phase 
and as an adsorbed species, we have optimized the struc- 
ture using a plane wave basis set. The optimized struc- 
ture is shown in Fig. [3] (a). The calculated W-C bond 
length of the free molecule obtained using a plane wave 
basis is found to be 2.06 A and the C-O bond length is 
1.16 A compared to 2.06 A^^ and 1.15 A respectively de- 
termined in experiment!^. Investigation of the electronic 
structure of W(C0)6 shows that the highest occupied 
state is dominated by W 5d orbitals. Indeed, the analysis 
of molecular orbitals also confirms this fact. The high- 
est occupied molecular orbital (HOMO) and the lowest 
unoccupied molecular orbital (LUMO) of W(C0)6 are 
shown in Fig. [31(b) and (c), which show that the HOMO 
is made up of W 5d while the LUMO is dominated by p 
orbitals of the CO ligands. 



C. Interaction of W(CO)6 with the fully 
hydroxylated Si02 substrate 

The FOH - Si02 surface corresponds to a substrate pre- 
pared under wet chemical conditions in the absence of an 
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FIG. 4: (Color online) Top panel: Three different configu- 
rations of W(C0)6 (CI, C2 and C3) considered for the ad- 
sorption on the FOH - Si02 surface are shown schematically. 
Similar configurations are considered on the POH - Si02 sur- 
faces where the CO ligands bond directly to the surface Si 
atoms. Bottom panel: Variation of the adsorption energy 
A^ for the different configurations with respect to the orien- 
tation of W(C0)6. Results are shown for the FOH - Si02 and 
POH - Si02 with and without inclusion of vdW corrections 



electron or ion beam. In accordance with the available 
experimental results, CO ligands of W(C0)6 have been 
used to bond to the substrate^^. Three different orien- 
tations are considered which differ in the number of CO 
ligands bonding to the surface. These configurations were 
subsequently relaxed without any symmetry constraints. 
A schematic representation of the three configurations 
considered (CI, C2 and C3; Cn involves bonding to n 
CO ligands) and their corresponding adsorption energies 
A^ are shown in Fig. |4](a) and (b). 

Our results illustrate that the W(C0)6 molecule 
bonded with two of its CO ligands to the substrate (C2) is 
the most stable configuration with an adsorption energy 
of -0.162 eV upon neglect of vdW interactions. However, 
the energy difference compared to the next stable config- 




FIG. 5: (Color online) Illustration of the role of the vdW 
correction in determining the bonding of W(C0)6 (C2) to 
the FOH - Si02 ((a) and (b)) and POH - Si02 ((c) and (d)) 
surfaces. The upper panel ((a) and (c)) correspond to calcu- 
lations without vdW corrections while the relaxed structures 
in the lower panel ((b) and (d)) include vdW corrections. 



uration, which involves three CO ligands bonding to the 
substrate (C3), is quite small (about 25 meV), which may 
be indicative for a certain capacity of the molecule to roll 
over the surface. In the relaxed structure, the minimal 
distances observed between the molecule CO groups and 
the surface OH groups are 2.07 A (CI), 2.14 A (C2) and 
2.06 A for configuration C3. 

The weak interaction between W(C0)6 and the FOH 
- Si02 surface results only in minor rearrangements of 
the substrate surface geometry. For instance, the Si- 
O bond distance on the adsorb ate- free surface is in the 
range of 1.62-1.65 A. For the relaxed configuration with 
adsorbate we do not observe any significant changes of 
these distances, illustrating the fact that the surface does 
not contribute significantly to the bonding. We also ob- 
serve only minor changes in the structure of the W(C0)6 
molecule bonded to the surface. In all the configurations 
we have considered, the W-C and the C-0 bonds that are 
involved in the bonding with the surface are shortened 
by 0.01-0.02 A which also confirms the weak adsorption. 

The role of vdW forces on these adsorption geometries 
is evaluated by adding vdW corrections to the calcula- 
tions which are not accounted for by DFT within the 
GGA functional. As expected, the absolute values of the 
adsorption energies are increased for all probed configu- 
rations (Fig. El (b)). For the most stable configuration 
(C2) we find an adsorption energy AE = —0.498 eV. 

We observe minor changes in the structural parame- 



ters when vdW corrections are included in the calcula- 
tions compared to the set of calculations without these 
corrections. For example, the minimal distance between 
the surface and the molecule in the three configurations 
are 2.11 A (CI), 2.20 A (C2) and 2.08 A (C3) and the 
changes in W-C and C-0 bond lengths are very small (of 
the order of 0.01 A) compared to the calculations without 
dispersion corrections. 

Although the minimal distance between the molecule 
and the surface are similar in both cases (C2 and C3), 
for the preferred configuration (C2), the remaining CO 
ligands are found to lie closer to the surface with the next 
nearest distances of 2.29 A and 2.65 A (see Fig. [s] (b)) 
compared to 2.82 A and 3.43 A (see Fig. [F] (a)) when 
the corrections are not implemented. This in fact has an 
impact on the orientation of the molecule on the surface. 
Furthermore, the differences between the adsorption en- 
ergies of the different configurations are now larger, thus 
indicating that vdW interactions tend to decrease the 
ability of the precursor molecule to spontaneously roll 
over the surface. 

The DOS of W(C0)6 on the FOH - SiOs surface for 
C2 is shown in Fig. |6] (middle panel) in comparison to the 
DOS of the slab without adsorbate (top panel). Results 
are shown for calculations including vdW corrections, 
which however have only minor infiuence on the elec- 
tronic structure (compare Fig. [5] (bottom panel)). The 
highest molecular level of the W(C0)6 adsorbate lies in 
the gap of the substrate, which pins the Fermi level close 
to the valence band. Thus the DOS of the Si02 substrate 
consisting mostly of O p states is shifted to lower energies 
with respect to the Fermi energy. Otherwise the DOS of 
the Si02 substrate shows only minor modifications with 
respect to the free substrate. The molecular levels of the 
W(C0)6 adsorbate are visible in sharp peaks and agree 
well with the energetic levels of the free W(C0)6 molecule 
obtained by our Turbomole calculations (also shown in 
the top panel of Fig.|6|. This further illustrates the weak 
interaction between the adsorbate molecule and the Si02 
substrate. 
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FIG. 6: (Color online) DOS of W(C0)6 on SiOs in com- 
parison to the clean Si02 surface. Shown are the total DOS 
(solid lines) along with the partial contributions of the adsor- 
bate (dashed lines) for the preferential configurations on the 
POH - Si02 surface (bottom panel) and on the FOH - Si02 
surface (middle panel). The top panel shows the total DOS of 
the Si02 surface (solid line) along with the energetic levels of 
the free W(C0)6 molecule. The Fermi energy is set to zero. 



D. Interaction of W(CO)6 with the partially 
hydroxylated Si02 substrate 

Partial de-hydroxylation of the surface of a wet chemi- 
cally prepared Si02 substrate is expected to occur under 
irradiation with an electron or ion beam or at elevated 
temperatures. To represent this situation in this study, 
we kept the Si atoms bonding to CO ligands untermi- 
nated, while all other dangling bonds of Si on the two 
surfaces of the slab were terminated by OH groups. The 
W(C0)6 molecule was placed in such a way that it bonds 
with one, two or three CO ligands to the unterminated 
Si surface atom and the substrate/precursor system was 
subsequently relaxed. 

Our results indicate that the W(C0)6 molecule is ad- 
sorbed more strongly on the POH - Si02 surface in all 



three orientations compared to the FOH - Si02 surface 
(see Fig. [4] (b)). Also on the POH - Si02 surface the 
precursor molecule prefers to be bonded through multi- 
ple CO ligands. In the multiple bonding configurations 
(C2 and C3) we observe a spontaneous fragmentation of 
the W(C0)6 precursor molecule with chemisorption of a 
W(C0)5 fragment on the surface upon removal of one of 
the CO ligands. 

The most stable configuration (C2) after relaxation 
with AE = —1.010 eV upon neglect of vdW interac- 
tions is shown in Fig. ^ (c). The fragmented CO lies 
at a distance of 4.15 A from the sub-carbonyl moiety 
(measured between 'W and CO) and 5.03 A from the 
surface (measured between C and nearest H atom of the 
surface OH group) and remains as a free CO molecule. 
The bond distances between Si and the O atoms of the 
CO ligands involved in the bonding are about 1.71 A and 



quite close to the Si-0 bond lengths of 1.62 - 1.65 A ob- 
served in various compounds^^, indicating a nearly op- 
timal interaction. The distances between the O of the 
two non-bonding CO ligands to the H atoms of the OH 
group are 1.94 and 1.98 A indicating the presence of a 
weak interaction which might also stabilize the molecule 
on the surface. 

A similar effect is also seen for the singly coordinated 
system (CI). When the W(C0)6 molecule is bonded by 
only one of its CO ligands to the POH - Si02 surface, 
the distance between Si and O of the nearest CO lig- 
and is 1.73 A with the fairly small adsorption energy 
AE = —0.218 eV. The distance of non-bonding CO lig- 
ands to the substrate are 2.31 and 2.58 A compared to 
2.08 - 2.48 A found for FOH - Si02 surfaces. In the three 
coordinated system (C3), where the W(C0)6 molecule is 
bonded by three of its CO ligands to the POH - Si02 sur- 
face, the Si-0 (of CO) distances range from 1.69 - 1.70 A 
with an adsorption energy AE of -0.723 eV. 

In the preferred configuration (C2) the W-C bond 
lengths of 1.87—1.94 A for the CO groups which are in- 
volved in the bonding to the substrate are shorter com- 
pared to the value found in the gas phase 2.05 A. On the 
other hand side, the bond lengths of the non-bonding CO 
ligands are slightly elongated (2.11 - 2.12 A), illustrating 
a weakening of these bonds. 

Calculations have also been carried out by including 
vdW corrections. No major changes are observed in the 
structural parameters of W(C0)6, (see also Fig. [s] (c) 
and (d)). However, differences are observed for the cal- 
culated AE on these POH - Si02 surfaces (Fig. [4] (b)). 
The most stable configuration (C2) has an adsorption 
energy A^ = —1.262 eV when vdW interactions are in- 
cluded. Moreover, the differences in AE between bi and 
tri -coordinate orientation of W(C0)6 on these surfaces 
becomes small, illustrating the presence of physisorption 
although dominated by chemisorption in these cases. Ir- 
respective of either the orientation or the inclusion of 
vdW corrections, fragmentation of W(C0)6 molecules is 
found to occur. 

The DOS of the preferential configuration of W(C0)6 
on the POH - Si02 surface (C2) is shown in the bottom 
panel of Fig. |6] Similar to the case of W(C0)6 on the 
FOH - Si02 surface, the overall shape of the partial DOS 
of the substrate is only little altered by the interaction 
with the precursor molecule, but shifted to lower energies 
with respect to the Fermi level since the highest occupied 
level of the adsorbate lies in the gap of the substrate. 
However, in contrast to the result for the FOH - Si02 
surface, the energetic levels of the adsorbate now have 
little in common with those of the free gas phase W(C0)6 
molecule, as a consequence of the stronger interaction 
and the fragmentation of the precursor molecule. 

In order to analyze the nature of the bonding between 
precursor and Si02 substrate in the adsorption process, 
we show in Fig. [7| the charge density of the highest occu- 
pied valence band of W(C0)6 adsorbed on FOH - Si02 
(see Fig. [7| (a)) and POH - Si02 (see Fig. [t] (b)) sur- 




FIG. 7: (Color online) Band decomposed charge density for 
the valence band maximum for the FOH - Si02 (a) and POH 
- Si02 (b) cases. 



faces. Our results indicate that the degeneracy of the 
W bd triply degenerate HOMO level is nearly preserved 
on the FOH - Si02 surface (cf. Figs. [3] (b) and[7| (a)), 
while it is strongly altered as a consequence of the sym- 
metry lowering while binding to the POH - Si02 sub- 
strate (cf. Figs.. [3] (b) and[7|(b)). This observation illus- 
trates the fact that the electron density on the tungsten 
is perturbed upon bonding to dehydroxylated Si atoms 
located on the Si02 surface which acts as an electron ac- 
cepting centre^^. Our computed Bader charges on the 
POH - Si02 surface before and after bonding illustrate 
that the dehydroxylated Si atom on the substrate gain 
0.9 electrons as a result of bonding. This interaction re- 
sults in the formation of a ligand-surface adduct which 
destabilises the molecule by weakening the W-CO bonds 
that were not involved in bonding to the substrate and 
ultimately eliminates one of its CO ligands. In fact, 
this behaviour is similar to the interaction of group VI 
metal carbonyls with the Al^+ sites of 7-AI2O3 result- 
ing in the formation of subcarbonyl entities bonded to 
the surface. ^^ Owing to the high adsorption energy, the 
formed sub-carbonyl moiety has less mobility on the sur- 
face than the unfragmented molecule. In summary, the 
observation of high adsorption energy of W(C0)5, the 
elongation of the W-CO bonds oriented towards the vac- 
uum and the fragmentation of one of the CO ligands on 
the POH - Si02 surface confirms that W(C0)6 bonds in a 
similar fashion in POH - Si02 the surface as on 7-AI2O3. 



IV. DISCUSSION: TOWARDS A REALISTIC 
SIMULATION OF THE GROWTH PROCESS 

The results obtained in the previous section show that 
W(C0)6 is only weakly adsorbed on a FOH - Si02 sur- 
face. The calculated adsorption energies in the range 
between AE ^ -0.1 to -0.5 eV may eventually allow the 
growth of W(C0)6 layers at low temperatures, but do 



not lead to the growth of metahic W deposits. Further- 
more, the smah energy differences between the probed 
configurations of W(C0)6 on the FOH - Si02 substrate 
are indicative of a high mobility on the surface, although 
the energy barriers between the configurations are hard 
to estimate. Experimentally, the formation of W(C0)6 
layers has been observed by exposure of a TiO2(110) sur- 
face to W(C0)6 molecules at T=140 K, which however 
desorb again after annealing to room temperature^. On 
the contrary, we observe chemisorption of W(C0)6 with 
a spontaneous fragmentation of the precursor molecule 
into W(C0)5 upon removal of a volatile CO ligand on a 
POH - Si02 surface. The possibility of mobility on the 
FOH - Si02 surfaces and fragmentation of the CO lig- 
ands on the POH - Si02 surfaces observed in this study 
illustrates that the molecule may roll over to the POH - 
Si02 surface regions and get activated for fragmentation. 
The rate limiting step in the decomposition of W(C0)6 
is believed to be the elimination of the first W— CO 
bond. In the gas phase, this elimination requires about 
45 kcal/mol^^ ^^, but is expected to be essentially smaller 
on the substrate and has been estimated to be about 
12 kcal/mol for W(C0)6 on Ni(lOO)^. Further, experi- 
ments have found that complete de-carbonylat ion o ccurs 
at around 150-205 °C in the CVD experiments^. This 
in fact, is supported by the weakening of bonds of the CO 
ligands not involved in the bonding to the substrate ob- 
served in our calculations, which simplifies the elimina- 
tion of further ligands by thermal or radiative processes. 
Our results thus show that the spontaneous fragmenta- 
tion of W(C0)6 on a POH - Si02 surface is a relevant 
process to be taken into account for the simulation of the 
initial stage of the growth of tungsten nano-deposits. 



- Si02 substrate results in a weak adsorption of the pre- 
cursor molecule. However, upon partial de-hydroxylation 
of the surface, we observe a chemisorption and fragmen- 
tation of the precursor molecule, which is explained in 
terms of the changes of the electronic properties and the 
charge transfer between the precursor molecule and the 
substrate. We observe that vdW corrections result in a 
re-orientation and stabilization of the molecules on the 
FOH - Si02 surface, but have only a minor impact on 
the adsorption geometry for the POH - Si02 surface. 

Our reported calculations provide insight into various 
aspects of the interaction of the precursor molecule with 
the Si02 substrate. In a real growth of nano-deposits as 
is the case of electron beam induced deposition of pre- 
cursor molecules on a substrate, there are a multitude of 
processes to be taken into account. The interaction of 
the beam with the precursor molecules in the gas phase 
for instance will lead to an activation or fragmentation of 
the precursor molecules. However, since in focused EBID 
experiments the fraction of gas phase molecules interact- 
ing directly with the beam is small, this process may be 
expected to be of minor relevance for EBID^^. Further- 
more, the interaction of the beam with the substrate will 
lead to heating and charging of the substrate as well as 
to emission of secondary electrons and photons. These 
processes can be partially taken into account by finite 
temperature molecular dynamics simulations, which we 
plan to do in the future but are beyond the scope of the 
present investigations. 
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